Does Good Mutation Help You Live Longer? 
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We shall argue below that our choice of an age- 
independent mortality leads to behavior which applies 
to systems with realistic mortality rates. The nature of 
our results also suggests that the details of the mutation- 
driven shift in offspring life expectancy is not crucial. 

Let C„(a,t) be the density of individuals with life ex- 
pectancy n > 1 and age a at time t. According to our 
model, the rate equation for Cn{a,t) is 
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C„(a,t). (1) 
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We study the dynamics of an age-structured population in which the life expectancy of an 
offspring may be mutated with respect to that of its parent. When advantageous mutation is 
favored, the average fitness of the population grows linearly with time t, while in the opposite case 
the average fitness is constant. For no mutational bias, the average fitness grows as t^'\ The 
average age of the population remains finite in all cases and paradoxically is a decreasing function 
of the overall population fitness. 

In this letter, we investigate the role of mutation on 
the age distribution and fitness of individuals within a 
simple age-structured population dynamics model. The 
basic feature of our model is that the life expectancy of 
an offspring, which measures its fitness, may be mutated 
with respect to that of its parent. While age-structured 
population models have been studied previously , rel- 
atively little is known about the role of mutation. When 
the individual reproduction rate is the fitness measure 
and population is regulated by externally imposed death, 
mutation leads to predominance by the fittest species [|j . 
In a related vein, it was recently shown that longevity is 
heritable |^,|| within the Penna bit-string model of ag- 
ing 1^. In these studies, the role of positive mutations 
was central. Our focus is quite different, as we study the 
dynamics of the fitness and age in a self-interacting pop- 
ulation as a function of the advantageous and deleterious 
mutation rates. 

When advantageous mutation is favored, that is, the 
offspring life expectancy (fitness) is greater than that of 
its parent, the fitness distribution of the population ap- 
proaches a Gaussian with average fitness growing linearly 
in time and dispersion increasing as t^/'^. Conversely, 
when deleterious mutation is more likely, there is a t~^/'^ 
approach to a steady fitness distribution. In the ab- 
sence of mutational bias, the fitness distribution again 
approaches a Gaussian, with average fitness growing as 
t^/^ and width growing as i^/"^. The average age of the 
population reaches a steady value in all cases and, sur- 
prisingly, is a decreasing function of the average fitness. 
Therefore within our model, a fitter population leads to 
a decreased individual lifetime. 

Our model is a simple population dynamics scenario 
which incorporates age structure and mutation. This dy- 
namics is based on the logistic model, N — bN — jN'^, 
in which a population with density N{t) evolves both by 
birth at rate b, and death at rate 7A'', with steady-state 
solution — b/j. The crucial new element in our 
model is that the life expectancy of each newborn may 
be mutated by ±t (with |t| = 1 without loss of general- 
ity) with respect to that of its parent. We also assume 
a constant age-independent mortality rate and birth rate 
for each individual. 

Each of these features represent idealizations of real- 
ity; for example, it would be more realistic to incorporate 
a mortality rate which is an increasing function of age 



The derivative with respect to a on the left hand side 
accounts for aging On the right hand side, the 

loss term jNCn accounts for death by competition and 
is assumed to be independent of an individual's age and 
fitness. As discussed above, the mortality rate is taken 
as age independent; the form C„/n guarantees that the 
life expectancy equals n. 

We account for the population of newborns as a bound- 
ary condition for C„(a = 0, t) [Q. An individual produces 
offspring with the same life expectancy at rate 6, and, due 
to mutation, produces offspring whose life expectancy is 
longer or shorter than its parent by ±1, with respective 
rates b±. Defining Pn{t) = daCn{a,t) as the density 
of individuals at time t of any age whose life expectancy 
equals n, then the boundary condition for C„(0, t) is 



C„(0, t) = bPJt) + b+P,,_i{t) + b_P„+i{t). 



(2) 



To determine the asymptotic behavior of the age and 
fitness distributions, it proves useful to first disregard the 
age structure and focus on fitness alone. From Eqs. (|^)- 
(0), the rate equations for P„(i) for n> 1 are 



P„+5+P„_i+6_F„+i, (3) 



with Po = 0. This describes a random- walk- like process 
in a one-dimensional fitness space which is augmented 
by birth and death due to the first term on the right- 
hand side. Using N{t) = I]P„(i), we find that the total 
population density obeys a generalized logistic equation 

dN , , ^ Pn 

= iB-jN)N-J2 — -b^Pu 4 



dt 



n=l 



where _B = + 6+ + 6_ is the total birth rate. 
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We now discuss the asymptotic behavior of these rate 
equations for three basic cases: subcritical - deleterious 
mutations favored (6_ > 6+); critical - no mutational 
bias (6+ — b-); and supercritical - advantageous muta- 
tions favored (6+ > In all three cases, the total pop- 
ulation density N and the average age A = N^^ J2^n, 
with An = J aCn{a) da, approach steady values. These 
are determined by a balance between the total birth rate 
B and the death rate 7-/V due to overcrowding. In the 
critical and supercritical cases, this leads to the steady 
state behaviors for the total density and the average age, 
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(5) 



The behavior in the subcritical case is more subtle, as we 
now discuss. 

Subcritical Case. Here a steady state is reached 
whose properties are found by setting = in Eq. (||). 
We solve this rate equation by introducing the generat- 
ing function F{x) = J2n>i -f" to transform the rate 
equation into the differential equation 



F' _ jN -b+1- 26_x 
'y ~ 6_ - (77V - b)x + b+x^ 



(6) 



Integrating Eq. (||) , subject to the obvious boundary con- 
dition F{1) = N, gives a family of solutions which are 
parameterized by the total population density N. To ex- 
tract a unique solution one has to invoke on additional 
arguments. First notice that N lies within a finite range. 
The upper limit is found from the steady-state version 
of Eq. §), {B - -fN)N = En-^Pn + 6-Pi > 0, to give 
jN < B. The lower limit is obtained from the physical 
requirement that all the P„'s are positive and therefore 
F{x) is an increasing function of x. From Eq. this 
leads to the inequality (7A^ — 6)^ > 46+6_. Thus 



be the basin of attraction for any initial condition with 
P„(0) decaying at least as fast as Conversely, an 

initial condition which decays as a" with a in the range 
X'^^^ = 1) should belong to the basin of attraction 
of the solution where, from the steady-state version of 
Eq. (^) in the large-n limit, the total population density 
is iV = {b+b-a+b+a~^) /j. We have verified this general 
classification of solutions numerically . 

Since the the steady state is approached exponentially 
in time for the classical logistic equation, TV = bN — jN^, 
one might anticipate a similar relaxation for our age- 
structured logistic equation (^. However, a numerical 
integration of the rate equations gives a power-law relax- 
ation of the total population density, Noo — N{t) ~ t^^^^, 
for a compact initial condition (Fig. |l|). This is also ver- 
ified by an asymptotic analysis of the rate equations 
A similar relaxation also occurs for the subpopulation 
densities with given fitness, P„(i). 
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b + 2y/b+b- < 7iV < B. 
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FIG. 1. N(t) versus t "^^^ in the subcritical case with b — 0, 
b+ = 1, — 2, and 7 = 0.5, and with the initial condition 
Pn(t = 0) = 0.1 for 1 < n < 10. The asymptotic inter- 
cept with the y axis gives the theoretically predicted value of 
iV,„i„ = 4^/2 « 5.6568. 



For any initial condition for the P„ with a finite support 
in n, only the minimal solution which satisfies the lower 
bound of Eq. (j^) is realized. This selection is reminiscent 
of the behavior in the Fisher-Kolmogorov equation and 
related reaction-diffusion systems 0|. 

To understand why the minimal solution is selected, 
consider the steady-state asymptotic behavior of P„ for 
n — > 00. In this limit, we may neglect the Pn/n 
term in Eq. (^. The resulting quasi- linear equation 
has the solution P„ = ^+^+ + ^-A", with X± = 

jN~b± y/{-fN - 6)2 - Ab+bJ\ /26_, and with A± < 

1. Thus the steady-state fitness distribution decays ex- 
ponentially with n. When the total population density 
attains the minimal value N-^i-^ = {b + 2-\/6+6_)/7, A+ 



achieves its minimum possible value A 
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where /i is the mutational bias. Since P„ 



the fitness distribution has the most rapid decay in n for 
the minimal solution. This minimal solution appears to 



For the relevant situation where the density N takes 
the minimal value, we integrate Eq. (^) to give the gen- 
erating function 



F{x) = N 



1 
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b+{fi - x){n - 1) 



(8) 



One can formally determine the P„ by expanding F{x) 
in a Taylor series. However, the asymptotic character- 
istics of the fitness distribution are more easily deter- 
mined directly from the generating function by using 



N Z^n=l — \ 
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ing this to Eq. (g), the first two moments of the fitness 
distribution are 
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The average age of the population may be obtained by 
first solving Eq. (|l|) in the steady state to give 
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(10) 



The average age then is ||] 
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(7iV)-i' 
x~ F{x) dx, 



(11) 



where the second line is obtained by using the expres- 
sion for C„(a) from Eq. (|l^) and the last line follows by 
expressing the sum in terms of an integral of the gener- 
ating function. The surprising feature that emerges by 
numerical evaluation of this integral (Fig. ^ is that the 
average age decreases as the population gets fitter! 
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FIG. 2. Average age A of the steady state population ver- 
sus /i^^ — y^b+ /6_ with 6 = 0, and with fixed total birth 
rate equal to one. For fi~^ > 1, the average age A = 1, while 
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Supercritical Case. When 6+ > 5_, the random 
walk in fitness space defined by Eq. (|^) is biased away 
from the origin and a continuum approach becomes ap- 
propriate in the long-time limit. Treating n as continuous 
and Taylor expanding the master equation for small de- 
viations about n, gives the following convection-diffusion 
equation, supplemented by birth/death terms, for the fit- 
ness distribution 
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(12) 



The difference between the advantageous and deleterious 
mutation rates now defines a bias velocity V = 6+ — 6_, 
and the average mutation rate plays the role of a dif- 
fusion constant D = {b+ + b-)/2. Integrating over all 
fitness values, the total population density obeys 



dN 



P{n,t) 



dn. 



(13) 



Since the fitness distribution is sharply peaked at (n) = 
Vt (see below), the integral on the right-hand side ap- 
proaches N/Vt. By setting ^ = in the resulting equa- 
tion, we conclude that 7A^ B—y^. This gives both the 
steady-state density, as well as the rate of convergence to 
the steady state. 

We now find the fitness distribution by substituting 
this asymptotics for N{t) into Eq. ( |l2|) to give 
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The birth/death term on the right hand side may be 
neglected, since (n) = Vt, and fluctuations about this 
average are of order t^/^ This approximation reduces 
Eq. (|lj) into the classical convection-diffusion equation, 
with solution 



P{n,t) = 



N 



exp 



{n)y 



ADt 



(15) 



This gives a localized fitness distribution with average 
fitness growing linearly in time, (n) = Vt, and width 
growing diffusively, a = ^2Dt. 

To determine the age characteristics, notice that 
asymptotically, the P„'s change slowly with time, so that 
the time variable t is slow. On the other hand, the age 
variable a is jast. Physically this reflects the fact that 
during the lifetime of a typical individual the change in 
the age characteristics of the population is small. Thus 
in the first approximation, we retain only the age deriva- 
tive in Eq. (|^). We also ignore the term C„/n, which 
is small near the peak of the asymptotic fitness distri- 
bution. Solving the resulting master equation and using 
the boundary condition of Eq. (0) we obtain 



C„(a,i)~P„(t)7iVe"^^'^ 
7iV2 



y/^TTDt 



exp 
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(16) 



Summing over the fitness variable, the total age dis- 
tribution C(a,t) = ^C„(a,t) is just a (stationary) 
Poisson, C{a,t) — •yN'^ e~'^^'^, and the average age is 
A = (7iV)~^ = B^^, in agreement with Eq. (||). 

Let us compare this average age to that in the sub- 
critical case; the latter is given by Eq. ( pi] ) with 7A^ = 
b + 2 1/6+6-. To provide a fair comparison (Fig. ||), take 
the total birth rate B to be the same in both cases. It can 
then be proved that the average age in the supercritical 
case is always smaller than that in the subcritical case . 
Individuals in a population with preferential deleterious 
mutations live longer than if advantageous mutations are 
favored! The continuous "rat-race" to increased fitness 
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in the supercritical case does not lead to an increase in 
the average life span. 

Critical Case. With no mutational bias, the fitness 
still grows indefinitely, but more slowly than in the su- 
percritical system. The equation of motion for P{n, t) 
is again given by Eq. (12), but with V set equal to zero 
and with N{t) is still described by Eq. ([l^). To derive the 
scaling behaviors of (n) and the width of the fitness distri- 
bution, we first use the fact that numerical integration of 
Eq. (|l^) again gives a localized fitness distribution. Thus 
we may estimate the integral on the right-hand side of 
Eq. (p|) as N/{n). This leads to 77V ^ B - Substi- 



tuting this into Eq. (12) for P{n,t) now yields 
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To determine the long-time behavior of this equation, 
we exploit the fact that the fitness distribution is peaked 
near n « (n). This suggests changing variables from 
{n,t) to the co-moving co-ordinates {y = n — {n),t). 
Eq. ( p7| ) then becomes 



dP 
'dt 



d{n) 
dt 



dP 
dy 



y 



P 



y 



P + D 



d^p 



(18) 



Let us first assume that the average fitness grows faster 
than diffusively, that is, (n) 3> Vi- With this assump- 
tion, the dominant terms in Eq. ( [T^ ) are 



d(n) dP 
dt dy 
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(19) 



These terms balance when {n)/{ty) ^ y/(n)'^. Using this 
scaling in Eq. ( |l8| ) and then balancing the remaining sub- 
dominant terms gives y ^ \/t. The combination of these 
results then give (n) - t2/3 rpj^jg justifies our initial 
assumption, (n) pg> \/t. Finally, writing (n) = (ut)^/'^, 
simplifies Eq. (|l^) to 
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(20) 



whose solution is the Gaussian of Eq. (|l_g), but with 
(n) = (itt)^/^. The value u = VSD is determined by 
substituting (n) = (ut)"^^^ in Eq. ( p^ and balancing the 
subdominant terms. To summarize, a Gaussian fitness 
distribution holds in both the critical and supercritical 
cases with the fitness distribution peaked at 
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supercritical case. 



(21) 



The age distribution in the critical case is obtained 
similarly to the supercritical case. The asymptotics of 



C„(a,i) is again given by a form similar to Eq. (|l^), 
which gives C{a, t) = ^N"^ ^-i^a g^f^^gj. gumming over n. 
Hence the average age is B"^, as in Eq. (||). 

While our discussion is based on a population dynam- 
ics with an age-independent mortality rate, this assump- 
tion does not substantially affect our main results. The 
crucial point is that old age is unattainable within our 
model. In the critical and supercritical cases, this is due 
to death by increased competition among fit individuals, 
while in the subcritical case, age is limited by the delete- 
rious mutational bias. Thus for a more realistic mortality 
rate which increases with age, similar fitness and age dy- 
namics to those outlined here would still result Q . 

In summary, in our population dynamics model, the 
average fitness grows linearly in time when advantageous 
mutations are more likely and the fitness approaches a 
steady value when deleterious mutations are favored. In 
spite of this fitness evolution, the average age of the pop- 
ulation always reaches a steady state. Intriguingly, this 
average age is a decreasing function of the average pop- 
ulation fitness. This paradoxical behavior arises because 
competition becomes keener as the population becomes 
fitter. 
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